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SOME INTERESTING FEATURES OF FREQUENCY CURVES 


By Ricumonp T. Zocu 


Introduction 


It is well known that in the normal error curve the points of inflection are 
equidistant from the mode. However it has never been pointed out that this is 
also a characteristic of all of the bell-shaped Pearson Frequency Curves. This 
fact can be most easily shown by placing the mode at the abscissa x = 0. 

Many rough checks have been developed for use in applying the Theory of 
Least Squares. The second part of this paper develops a rough check on the 
computation for use when fitting a Pearson Frequency Curve to a set of observa- 
tions. No rough checks on computation are given in textbooks on Pearson’s 
Frequency Curves. 

At present it is customary to follow a separate procedure for each Type of 
curve when computing the constants of a Pearson Frequency Curve. The 
third part of this paper shows how a single system may be followed for all Types. 
A single procedure is very desirable in order that the rough check of Part 2 may 
be quickly applied. 


Part 1. Points of Inflection 


Perhaps nothing brings out the limitations of the beli-shaped Pearson Curves 
in a more striking manner than a discussion of their points of inflection. In 
dealing with frequency curves it is well known that any curve can be fitted to a 
given distribution and that the real problem in curve fitting is the selection of a 
curve. Figures 1, 2, and 3 illustrate three hypothetical histograms. All three 
of these histograms are bell-shaped yet none of them will be closely fitted by 
any of the Pearson Curves. The reasons will be pointed out presently. 

The differential equation from which Pearson derived his system of frequency 
curves is : 

dy y(x — P) 


dx bax? + bx + bo 


By putting  — P = X, i.e. by placing the mode at the abscissa X = 0, this 
differential equation may be written: 


dy _ yX 
dX +BX+BX+ By 


where the + or — sign is taken according to the type of the curve. (It will be 
shown later that the constant term of the denominator must be less than zero.) 
1 
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Since in the Type III curve B, is 0 and in the ““Normal Curve’”’ both Bz and B, 
are 0 it will be advantageous to consider the general case of 

dy yX 

dX F(X)’ 
where F(X) is an integral rational function of the nt degree, at once rather than 
considering special cases first. 


If 
dy yX 


dX F(X)’ 
then 
d*y y 


7X0 = (FORE + PO) ~ XPOO). 


9 
« 
. 


In order to locate the points of inflection, Xi is equated to zero. Then we have: 


X? + F(X) — XF’(X) = 0. (1) 


This equation is always of the same degree as F(X) except when F(X) is linear or 
constant. Hence we have proved the Theorem: If y = G(X) be the solution 
of the differential equation 


dy yX 


dX F(X)’ 


then the number of points of inflection of y cannot exceed the degree of F(X) 
when F(X) is of degree greater than one. 

Now F(X) = B,X" + B,iX"*" 4+ .-- + BX? + BX + Bo. Whence 
equation (1) can be written in the form: 


(1 — n)B,X* + (2—n)B,-1.X*" + (3 — n)B,-2X"? 4+ --- 
+ (r —n) Bary X**+1 4 --- — 3B X4 — 2B,X* + (1 — B.) X?4+ B= 0. 
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Hence we have established the Theorem: The coefficient of the linear term of 
X in the equation of the points of inflection is zero. 


mu 


Fie. 3 


For the ‘‘Normal Curve ” and also for Type III, 
B, = B; = By = tee = B, = 0. 


Hence the points of inflection of these two Types are given by X = ++/ — Bo. 
For Types I and II, B is positive and Bz; = By = --- = B, = 0, and the 
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points of inflection are X — Hence the points of inflection are 
— 2 
undefined if Bz = 1, are pure imaginary if B, > 1, and real if B, < 1. 
For Types IV, V, VI and VII, Beis — and B; = --- = B, = 0, and the 


oints of inflection are at X = wal ‘ 
‘ i+ “TB 
In some of these Types it may happen that the abscissae of the points of 
inflection though real will lie beyond the range of the curve. Thus Types III 
and VI seated — 1 or 2 points of inflection, the single point of inflection occur- 








ring when | /= + B, > the range of the curve in the direction that the range is 
2\ 


limited. Type if, may have 0 or 2 points of inflection, as there will be no real 
points of inflection when B, 2 1. Type I may have 0, 1 or 2 points of inflection. 
Types IV, V and VII as well as the ‘‘Normal Curve’’ always have 2 and only 
two points of inflection. 

Now it should be noted that when one of the eight bell-shaped Pearson curves 
has two points of inflection then the abscissae of these 2 points of inflection are 
equidistant from the abscissa of the mode. In figure 1 a point of inflection will 
be at abscissa b and another at abscissa a. (M is the abscissa of the mode.) 
Since b — M ¥ M — anone of the Pearson curves will fit this histogram closely. 
In figure 2, points of inflection occur at abscissae a, b, and c. Since a Pearson 
curve can have at most two points of inflection no Pearson curve will fit this 
histogram closely. In figure 3 there are four points of inflection and no Pearson 
curve will fit this histogram closely. 


Part 2. Range 


DEFINITION: A beli-shaped curve is a continuous curve which starts at zero 
(or zero as a limit), rises to a single maximum, at which maximum point the 
first derivative is zero, and then falls to zero (or zero as a iimit). 

Or, more formally, y = G(2) is a bell-shaped curve if G(z1) = G(x2) = 0 and 
if G’(P) = Oand G’’(P) < 0 where G(z) is continuous and does not vanish in the 
interval from 2, to x2 and P is a unique point in this interval. 

If a bell-shaped curve has the value of zero at two finite points, one on each 
side of the maximum (mode), it is said to be of limited range in both directions, 
or briefly, of limited range. 

If a bell-shaped curve has the value of zero at only one finite point it 1s said 
to be of limited range in one direction, or also of unlimited range in one direction. 

If a bell-shaped curve has the value of zero only at + ©, i.e. at no finite points, 
it is said to be of unlimited range in both directions, or briefly, of unlimited range. 

THEOREM I: If F(x) can be separated into a finite number of factors each 
either of the form (x — r,) or (x? + 2r; x + rj + 19,) where no real root is 
repeated and y = G(z) is a bell-shaped curve which is a solution of the differential 
equation 
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dy _ y(t — P) 


dx F(z) ’ 
then if F(x) has no real roots, y is of unlimited range in both directions; if all of 
the real roots of F(x) lie on the same side of P, y is of limited range in one (that) 
direction; if at least one real root of F(x) lies on one side of P and at least one on 
the other side, y is of limited range in both directions. 
Proor: If F(x) = 0 when x = P, we have 


dy _ y 
dx g(x) 


where g(x) = F(x) + (« — P). This derivative is zero only when y = 0 or 
g(x) = +. Hence the solution does not have a finite maximum and therefore 
is not a bell-shaped curve. If F(x) > 0 when z = P, we have 


d?y , 
fy. ph | @— Pet Fe) -— @- PS Pe | 


d?y 
dx? | ,=p = Fy 


which is greater than zero ae since at a maximum the second derivative must 
not be greater than zero, in this case the solution would have a minimum at 
x = P and therefore would not be a bell-shaped curve. As the theorem concerns 
only those solutions which are bell-shaped curves, F(x) < 0 when « = P. If 


[F(z)] 


; d ; 
F(x) = 0 when xz + P then oe = + ~ unless y is also zero. Assume y + 0. 


Since F(z) is negative, if y ~ 0 when F(x) = 0 then os — — » as F(r) > 0, 
for an x > P, and changes to + © as F(x) changes sign on passing through the 
value 0. Hence the curve would contain another maximum before falling to 
zero and therefore the solution is not a bell-shaped curve. Similar reasoning 
holds for an x < P. Therefore if y + 0 when F(x) = 0, the curve is not bell- 
shaped. If y = 0 when F(x) = 0, the curve has its range limited at this point. 
That is, any real number which makes F(x) vanish will also make y vanish if y 
represents a bell-shaped curye. Hence if all of the real roots lie on the same side 
of P the curve is of limited range in that direction only, while if at least one of 
the real roots lies on each side of P the curve is of limited range in both direc- 
tions. If F(x) contains no real roots it does not vanish for any real value of z. 
In this case, by partial fractions the differential equation becomes: 

dy _ k,, dx ky, dx Qk (x + 1,) dx 


y («+r 4ri, (+n)? +75, (x+r)?4+ 7%, 


Qko(x + r,) dx 
(2+r,)? +76, 
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On integrating, 


zt+r, 
ki. are tan 


y=C [ (@ +r)? + wi [ (« + 1)? + ro, |‘ ine il "Oy 





Hence y does not vanish for a finite real value of z and the Theorem is fully 
established. 

TueroreM II: If F(x) can be separated into a finite number of factors each 
either of the from (x — rj) or (x? + 2rja +r} + ro;) where no real root is repeated 
and y = G(z) is a bell-shaped curve which is a solution of the differential equation 
¥ = ey? , then if y is of unlimited range, F(x) contains no real roots; if y 
is of limited range in one direction, all of the real roots of F(z) lie on the same 
(that) side of P; if y is of limited range in both directions, at least one of the real 
roots of F(x) lies on one side of P and at least one on the other. 

Proor: By partial fractions the differential equation may be written: 








a a re 

y t—Tly t— N32 (x + 1)? + TO, 
eda, hee + ry) dx | Bhgo(e + To3) dx 
(x + 19)? + rj (x + 3)” > TO, (x + To)” + Ts 





and on integrating: 





Tre 















k ° 2 k kay are tan — 

, “19 © 2 8 r 

y = C(x — ry) (a — ry)A2 --- [ (x + 1)" + ro, | m ce @ 0} 

Hence y = 0 for x = ru, 72, --- and for no other finite values of x provided ky, 
ky, --- are positive. If one or more of the k;; are negative, y = © at such 


points and unless some r;; closer to P has previously made y vanish, the curve 
is not bell-shaped. Therefore, for bell-shaped curves, the exponent of the factor 
containing the real root of smallest absolute value on each side of P is positive. 
Therefore: if y is of limited range in both directions, at least one real root lies on 
each side of P; if y is of unlimited range in one direction, all of the real roots lie 
on the same side of P; if y is of unlimited range it contains no real roots. Hence 
the Theorem is established. 

The effect of repeated real roots will now be considered. If a real root is 
repeated an odd number of times at x = r, then F(x) changes sign at x = r 
and the first theorem is true. Ifa real root is repeated an even number of times 
at x = r, then F(x) does not change sign at x = r and we know that either (a) 






y = Oatz =7;or (b) yis finite and + Oand = +o atx =7, i.e. there is a 


point of inflection at x = r. It will now be shown that (b) cannot occur. If 





, . : d ‘ 
case (b) is possible, y is continuous at x = r, 7 = + « according as (r — P) $ 0 
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d — s . 
moreover c does not change sign in the neighborhood of the point x = r, and 
7 
dy 4 : s < 
dx changes sign from + x to — © or vice versa according as (r — P) §$ 0. 
dx? 


Now 


d*y y 2 d 
7,3 = taal — P)’ + F(x) — (x — P) dz Fa) |. 


ay 
dx? 
possible to select a neighborhood such that 


Whence if y is finite and + 0, does not change sign at x = r because it is 


(2 — P| >| Fe) — @ - P) Fe) 
| ar 


for an x differing from r by e where « is a small positive quantity. Therefore 
case (b) is not possible and y = 0 when a real root is repeated an even number of 
times. That is to say the range of the curve is limited at a point where a real 
root is repeated an even number of times. Thus Theorem I always holds for 
repeated roots. 

For Theorem II it is clear that this Theorem holds for repeated roots when a 
non-repeated root lies closer to P, and on the same side, than the repeated root. 
Suppose that the repeated root is the nearest root to P (on a given side of P). 
Then by partial fractions: 


dy_ _kude , kydr | kudr , 4, kadr | kedt 
(a = ra) (x — 42) 


y (t—rn) | (f—rn)? | (x — ru)? 


ig 50 a ape ee kde, 4 Phe + tn) de 
(z+) +70, (2 + rao)” + To (x + 7)? +75, 


and on integrating: 


y = C(z — r,,)"(z — ry)*"(z — ree)* --- [ (@ + 1)* + ro, |"* 


2t+r k k 
ko, are tan Se pe xia Da es 


+ — 
- "Oy (2-rn) = 2(2—-r 4)? 


Hence y can = 0 only for z = ry or for x = ra, Te, --- and for no other finite 
values of x. Since by hypothesis y is bell-shaped, then the proper k;; must be 
positive and Theorem II always holds for repeated roots. 
Theorems I and II can now be combined and generalized in the form: 
TuHeEoreM: If F(x) is a polynomial with real coefficients and y = G(z) is a 
bell-shaped curve which is a solution of the differential equation 


dy _ y(x — P) 
dx F(x)’ 
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then the necessary and sufficient condition: that y be of unlimited range in both 
directions is that F(x) have no real roots; that y be of limited range in one direc- 
tion is that all of the real roots of F(x) lie on the same side of P; that y be of 
limited range in both directions is that at least one real root of F(x) lie on one 
side of P and one on the other. 

Coro.uary: F(x) must be negative throughout the range of y. 

Suppose now that we have some statistics which we wish to graduate and the 
statistics are of such nature that we would expect a bell-shaped curve, rather 
than a J- or U-shaped curve, and we desire the best fit: If we use a curve which 
is a solution of the differential equation 
dy _ y(x — P) 


dx —S- F(x) ~~ 
(the Pearson Curves being special cases) to fit the statistics and if in computing 
the constants for the curve one of the following cases arise: 
(a) by < 0 when this constant is computed, 
or (b) Bo < 0 when the origin is moved to the mode, 
or (c) arootis located within the range of the statistics then it means that: 

1. A mistake may have been made in the computation: thus the Theorem 
just established provides a rough check on the work of computation, 

2. If no mistake has been made in the computation it may indicate that the 
bell-shaped Pearson Curves will not closely fit the statistics and that some 
other graduation curves be used, e.g. the Gram-Charlier Types A or B might be 
tried, 

3. If no mistake has been made in the computation it may happen that one 
of the bell-shaped Pearson Curves will give an excellent fit but a different method 
than or a modification of the Method or Moments should be used in order to 
compute the constants. 





Part 3. Computing the Constants 


At present, the constants of a frequency curve are computed as follows: 
First the moments are computed about an arbitrary origin, then the moments 
about the A.M. are determined, then 6; and 2 and the criterion are computed, 
after which the type of curve can be selected. From this point a separate 
procedure is followed for each curve. Now in the above method one will not 
know whether a root has been located in the range of statistics or not. 

Take Pearson’s differential equation 


dy ye P)_ 
dx box? + yr + bo 
Put X =x—P. ThendX = dzandz = X + P, and 


dy _ yX a yX 
dx bo(X+ P)?4+0(X +P)+ bo  b.X2+ 2PWX + bX + P%. + Pb, + bo 
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nh 
2Pbe + bi = By 
P%, + Pb; + bo = Bo. 


Then we have 


dy _ yX - dy y(a— P) " (1) 
dX B,X?+ BX + By dx B(x — P)? + Bia — P) + Bo 


It should be noted that for a particular curve, Bz, B; and Bo are constants; 
i.e., their values do not change with a change of the origin. The values of }; 
and by do change with a change in the origin. 


If we clear equation (1) of fractions, multiply by e’ and integrate with respect 
to x over the range from 2; to x2, where 


Aon2 han? 
meta tate d 
€ 2! = e”* ydz , 
z 


1 


then successively differentiate with respect to », and equate coefficients of 
like powers of n, we finally obtain: 


1 — P+ B, — 2PB, + 2B, = 0, 

ho + Bo — PB, + P?B. + Bix — 2PBr, + 3Bodo + Bord? = 0, 
As + 2d2B, — 4PBods + 4Bods + 4Bdsrho = 0, 

hs + 3BiA3 — 6PBodrA3 + 5Body + 6BodAzZ + 6Brdds = 0. 


Since we can compute the moments from the raw statistics and the semi- 
invariants from the moments, we may regard ie, A3 and A, in these equations as 
knowns and the Bo, B;, B2, P and \; as unknowns. But the origin has not yet 
been specified. Let the origin be placed at the A.M. where uw, = 4; = 0. As 
Xo, Az, Aa, Bo, B; and Be are unchanged by a change of origin, we have: 


B, — Py — 2PoB, = 0. 

he + Bo — PoB: + P2B. + 3Bad. = 0, 

As + 2Bire — 4PoBodo + 4Bod3 = 0, 

hs + 3Bid3 — 6PoBodr3s + 5Bodrs + 6BodAzZ = 0. 


by = Bo — PoBi + PiBzo, | 
b; B, _— 2PoB2, 
b, Bs ; 
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then 


b —P=0, 

de + bo + 3b,r2 = O, 

As + 2b re + 4b,A3 = 0, 

hs + 3b; As + 5dr + 6,432 = 0. | 


By reversing the transformation (4) we get: 


B, = by 


- 


? 


/ 


B, = b; + 2Pob, (6) 
By = by + Po(bt + Pobs). J 


Now the above theory suggests the following procedure for computing the 
constants of a frequency curve: First the moments are computed about an 
arbitrary origin, then the semi-invariants are computed (or alternatively the 
moments about the A.M., either step involves about the same amount of work), 
then the equations (5) are solved and then by means of equations (6) the Bs, 
B, and By are computed. Next solve the quadratic equation 


BX + BX + B = 0. 


The character of the roots of this equation indicates which type to use and it is 
unnecessary to compute the criterion. The constants of the frequency curve 
are simple functions of the roots of the above quadratic equation and can be 
readily found by integrating the diff. eq. (1) being careful to write the solution 
as a function of X = x — P. The rough checks mentioned in Part 2 can be 
quickly and conveniently applied when this procedure is followed. 


JEORGE WASHINGTON UNIVERSITY. 





A RECONSIDERATION OF SHEPPARD’S CORRECTIONS 
By W. T. Lewis! 


In computing the moments of a frequency distribution it is customary to find 
first what are known as the raw moments. These are obtained on the assump- 
tion that all the material of each class interval is concentrated at the middle 
point of the interval. It introduces what is called a grouping error because in 
fact the material does not all lie at the middle point. To compensate for this 
error W. F. Sheppard? derived a set of corrections. The hypothesis underlying 
his method is that the distribution may be regarded as similar to one to which 
the Euler-MacLaurin summation formula without its end terms may be applied. 
He presupposed such a curve, found its true moments, and then the raw moments 
that would be obtained if its area were concentrated at several equidistant 
abscissae. The relationship between these raw moments and the true moments 
of the curve furnished him with the corrections required for that distribution. 
If now our observed distribution may be supposed to be sufficiently like that one, 
we may use his corrections also on the observed data. One may note four points 
of criticism. 

(1) The given distribution may not be similar to the one suggested, in the 
sense that it would be close to such a curve if the intervals of grouping were 
made very small; or at all events the purpose of finding the moments may be in 
part to decide whether or not it would become such a curve, and so one would 
not like to assume that to be true at the outset. A special case of importance 
in which this last is true occurs when one is finding the moments of a sample in 
order to determine whether it may have been drawn from a presupposed universe. 
It is inexact to use raw moments but it is illogical to use corrections that have 
been proved only for the universe being tested. 

(2) Sheppard’s argument does not make use of the one certain fact that is 
given in the hypothesis, viz: that the partial area of the given distribution over 
each class interval is exactly as stated. In fact, if, following the argument of 
some authors, the given curve be assumed to be exponential, it obviously cannot 
have partial areas everywhere exactly equal to the several given frequencies, 
for in particular its partial area is not zero beyond the given range. 

(3) It is common to find distributions which do not have high contact at the 
ends of the range and for them Sheppard’s corrections certainly fail. To 
obviate this criticism new corrections have been derived by Pairman and Pear- 


1 With the assistance of Burton H. Camp. 


2 The true values are given on page 220 of ‘‘Mathematical Part of Elementary Statistics, 
by Camp, D. C. Heath and Company, 1931. 
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son for the so-called abrupt cases. These new corrections are adequate to care 
for the abrupt cases but involve so much computation that it is a fair question 
whether it would not be simpler, first to distribute the given material over each 
interval by a smoothing process, and then to find without corrections the 
moments of the smoothed distribution. 

(4) Even if one admits Sheppard’s method in general, waiving the dubious 
question as to whether it is proper to start with an assumed curve instead of 
starting with the given distribution, it is doubtful whether there are any curves 
which have exactly the properties required. The high contact hypothesis may 
be put in different language as follows: using the notation of the Handbook* 
page 92, let f(x) be the curve and x, be the middle point of the slice. It is 
assumed that 


—o 


> xis) = [ x’ f(x) dz; 1=0,1,---; r=0O0,1,---; 
t 


c being the class interval. This means that if the moments of the curve be 
found by using mid-ordinates times class interval, instead of areas, one will obtain 
exactly the true moments of the curve, and that this will remain true for all the 
curves which are derivatives of this curve. This property is certainly not true 
of the normal curve; but it is almost true when 7 and the class interval are both 
small, and it is probably due to this fact that Sheppard’s corrections seem to be 
good in practice. 

Moreover, this high contact hypothesis cannot be true for any function over a 
limited range if the function is developable in Taylor’s series about one end of the 
range. For the only function which has the required properties is identically 
zero, since the function and all its derivatives are required to vanish at that end 
of the range. 

The primary purpose of this paper, therefore is to derive corrections similar to 
Sheppard’s with a different set of assumptions. The results may be used as an 
approximate substitute for both Sheppard’s and Pairman’s. That is, they will 
apply approximately to both extreme cases and to the intermediate cases; on the 
whole they give better results than Sheppard’s and are not so difficult to admin- 
ister as Pairman’s. , 

The argument runs as follows. When a distribution is given merely by class 
intervals, there is no way of knowing exactly what the distribution would have 
been had the class intervals been smaller; we do not know that we have a sample 
from an exponential curve, and even if we did we would not know that this 
sample would lie close to the exponential in form. We shall, however, try to 
draw a graduating curve in such a manner that (a) its partial area over each class 
interval will equal the frequency of the given distribution over that interval; 
and (b) its form within each class interval will be such that it will pass smoothly 
into the adjacent portions to the right and left. A good way to do this is by a 






2H. L. Rietz, ‘Handbook of Math. Stat.”” Houghton Mifflin Co. (1924). 
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freehand graph, frankly recognizing that there are many forms that will do 
equally well. To obtain a numerical result it is necessary to use the equation 
of some curve. Again frankly recognizing that there are many types which 
will do equally well we choose the simplest to handle: 


y=a- bt + ct?. 


Let the relative frequency distribution be defined by f(z), —m Si Sn, m,n,t 
being integers. To satisfy (a) we have the equation 


[ova = fi). 


To satisfy (b) we shall let 

y= Af@) +f/@4+ lift =7+3. 
The latter will hold for all values of 7 from —m to n — 1 inclusive, but the end 
intervals require special treatment. Here in order to satisfy as well as possible 
both the high contact and the abrupt cases, we wish to let the material be 
distributed according to the way the curve is behaving over the two nearest 
intervals on the right (at n) or left (at —m) rather than by the addition of zero 


frequencies beyond the given limits. To do this we let the slope of the para- 
bolas be zero at the extremes: 


dy 9 


ae att = —m—jandt=n-+}. 
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Then, if for example the frequencies are increasing as one nears the right end 
interval, the curve will rise over the right end interval; if they are decreasing, 
it will fall. These three conditions are sufficient to determine a continuous 
curve of the sort indicated in the figure. The exact moments of the curve may 
be found by integration and expressed in terms of the raw moments. The 
details are tedious and of an elementary nature and will be given only for the 
mean value 7;. 

To determine the coefficients of the parabola y = a + bt + ct? for the rectangle 
at t = 7 we may write the following three equations; the first complying with the 
requirement that the area under the parabola from ¢ = 7 — } tot = 7 + 3 equals 
the area of the rectangle at ¢ = 7, the second and third giving the ordinates at 
@— Zandi + 3 respectively: 


f@ = fe (a + bt + ct’)dt, 


LOAIETD aed DteG4+D, 


te SS 0e8t = Bae ~e. 














Solving these three simultaneous equations we get for a, b, and ec: 


=~ ¢- a0 + (8 -5- cans (B+i-Dav-o, 


g 


b = Gif) + G — 3)f@ 4+ 1) — G4 38)f@- 1D), 
c= — 3X) +3f/G+0)+2f@—-1), 


and these hold for —m +1575 n-—1. 


For the parabola y = a, + bjt + ct? over the first rectangle, i.e., where 
1 = —™m, we get the equations: 






f(-— m) = i. (qj + bt +c dt, 


m—%Fs 





f(— m) — m + 1) = a, +bi(— m+ 4) +c, (— m+ 3)?, 








bi + 2c, (— m — 3) = 0, 


and their solutions: 





a = ¢ (m? + m — ys) f(— m+ 1) — 3 (m+ m— 35) f(— m), 
by = § (2m + 1) f(— m+ 1) — 3 (2m + 1) f(—m), 
a = 3f(— m+ 1) — 2f(—m). 






= 
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Similarly for the parabola y = a, + b,t + c,t? through the last rectangle at 
1= nwe get 


f(n) = fo (a, + bat + c,t?)dt , 


fin) + fn — 1) 
2 


b,+2en+e=0, 


= a, + b, (n — 3) +, (n — 9), 


and for the constants 
= i (n? +n — ye) fn — 1) — §(? +n — 35) f(n), 
= — $(1 + 2n) f(m — 1) + ¢ (1 + 2n) f(n), 
ifm — 1) — fm). 


Having obtained the constants for the graduating curve we will determine 
the moments of this curve in terms of those of the given frequency distribution. 


Notation: Let the class interval be c = 1; let », = > i f(i) be the uncor- 


rected st* moment of the given frequency distribution about the given origin; 


let ws = = (¢ — »)*f(z) be the uncorrected st* moment of the given fre- 


quency distribution about its uncorrected mean; let 7, be the corrected value of 
the st* moment about the given origin; and let 7, be the corrected value of 
the st moment about the corrected mean. Thus », and yu, apply to the rec- 
tangles, and 9, and ji, apply to the curves as follows: 


=-~% 


2, J 


i+} —m+} 


(a + bt + ct?)dt + | t*(a; + bit + c,t?)dt 


i, = 
——e | 

“ ™ 
1 3 


—i 
2 


n+4 
+ [ t*(a, + b,t + c,t?)dt, 


Jn—} 


n—1 
—m+} 


Ry. sm Zz [- ; (t — p,)* (a + bt a ct?)dt Sal (t _ p,)* (a, + bit + ct?) dt 


—_—ee | 
m 
i=—m+1 2 


n+} 
+ (i = D,)* (a, + bat + c,t?)dt . 
n—4 


Using these symbols we have for the first moment about the given origin: 


ne 


ta, + byt Leyt?)at 


. i+} 
Zz i t(a + bt + cf?)dt + 
t=—m+1 ; 


2 as 


fn+h 
+ | . t(a, + bat + cpl?)dt 
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oz 


- Seso(oa aed) 

+ | - a,m-+ by (m: + 5) — € (ms +- ™) | 

+ | n+b (x 5) c (n+ ")| 
nN + n + 12 + n 4 . 


Substituting the values for the constants this becomes 





ent = ‘i Ie _ 32) fli) + (= “ae :) fi +) 
37? 2 1 , 
‘ & “ie 5) i 7 » | 
+ (2 + a) 6@) + — 3064) — 4430/6 — DI 
¥ (* + ‘) [— 3) + 8G + 4356 — vi} 


+ {— m[¥ (m? + m — Ps) f(— m+ 1) — 2 (m? + m — 15) f(— m)] 
+ (m? + yz) [F (2m + 1) f(— m + 1) — § (2m + 1) f(— m)] 


al (m* + *) [3 7(— m +1) — 3f(— m))j 
+ {nlf n? +n — ah) fln — 1) — 2 42 — 3D fad) 
++) | - 3 (1 4 an) sin -1) + 3 + 2n) s(n) | 


s(esi)[ie—0 ~ th 


a3 


= a [io +3 + 5 fE+1) — 5 SG - | gi(- m+) 


— (m+ ii)A- ee 5, In —1)+ (" + ig) 


Zz f(i) = » if(t) — (— m) f(— m) — nf(n) = » + mf(— m) — nf(n). 


—m+1 —m 


sh z fi +N = 3 > fk) = xe - aim +1) — f(—m) 









1 1 1 
*_-gh- a+) - 2 K- m) . 
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a1 n—2 


a z fi — ~ 94 7] f(j) = aD fj) - gin _ — 544) 


i=—m+1 


1 1 
= 54 fm —1)- aI) : 
1 ] 1 1 
= vy, + mf(— m) — nf(n) + —- 54 I\- m+ 1) — at m) — 54 


+ a f(n — 1) + 54 1") + 5 Al- m + 1) 


<n (m 4. - 4) m) — ag ft — 1) + (n +. X) f(n). 


5 
hen-oiK- m) + 5 * fn) + 5 5 m+ 1) — pim—0). 
Using this same notation and method for the higher moments we get 
1 _9 5m 7 5n 7 
ae + (= . z)K- m) + (= + x) f(n) 
+($- 1) S- m + 1) +(3 1) i 1) 
24 ~ 240 +t 24 ~ 240)°°" 

ie ma ee 
= v3 — 3Pifle 7% + f(— m) lz m? — om — = 


+40] 3 n +5 = “n + in| + f(— m+ 1) Ek + 5, + aa | 


—n 7 1 
—* ost -.ae 


= 14 — 4737; — 6f205 


5m 2 Pat 17n 313 
+ f—m)| m4 ) 5; 40 * 30 ae | 


eftem eof 52 


SPECIAL CASES 
The above formulae are rather long and in practice the special cases below 
will frequently be preferred. 


(a) We may usually take the origin at or very near the middle of the range so 
that m = n, at least approximately. 





18 
If m = n: 
Zim) + 2 $n) + FeS(—m +) — Bsn — 2). 


ji2 = v2 — : — i+ (= + a )l(—m) + f(n)] 


+ (SP Ap)i—m + 4 se — 


f= ve — Brita — 7% — oh + | Mg Bm a IT Uegeny — s(—m)] 


+|% + 5 m+ aay [—m +» —f(n —1)). 


ae -~ 
4 He Py 17 


fy = mu — Aisi — 69,5; — 1 -F—- 5 - we 





—m mm m 
+| 12 — 7 ~ 307 aa li- m+1)+f(n — 1)] 





| = 21m? 17m 313 


12 40 30 + eb (—m) + f(n)]. 


(b) Except in the abrupt cases the end frequencies and the difference between 
those next to the ends will be so small (relative to unity) that they will have a 
negligible effect on the corrections. If m = n as in (a), and if also 


f(—m) = f(n) = 0 and f(—m + 1) — f(n — 1) = 0: 


y= WN. 


i 1 —m 1 
a + f(- m+ 0 - al 


35, i V1 =$ 
wae SS — es 


These formulae have been written in the form which makes the computing 
simple. The following makes a comparison with Sheppard’s corrections easy. 
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Me — a3 + f(—m + Y 52 
= ws + (2 + a) s—m +1). 


Me 
= — 8 Bs mg a 


The following special case is also useful in comparing my formulae with 
Sheppard’s. 

(c) Let f(—m) = 3 f(—m + 1) and f(nm) = 3 f(2 — 1). This produces a 
graduating curve which is exactly tangent to the t-axis at the ends of the range 
and is everywhere continuous—though it does not have continuous derivatives 
at certain isolated points. It is, however, a curve which to the eye cannot be 
distinguished from the type assumed in the Euler-MacLaurin theorem, which 
lies at the base of Sheppard’s formulae. My corrections become: 


y= VN, 


ee a + 5 (—m) + f(n)}, 


*i—m) + son + |" — 2 |p—m + [244 |sem, 


_ 7~ a + | -- m)? +nut+m+= = | m) 


+2| —n?—-n+n+ oe [seo 


Sheppard’s are: 


is= ay —2 t+ oo 


Let us compare my results with Sheppard’s in the very special case in which 
f(—m) = f(n) = 1/7, f0) = 5/7, m = n = 1. The odd moments vanish. 
My corrections for ue and ps are 


i, = 0.2214, jy = 0.1870. 
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Sheppard’s are 
fe = 0.2024, fs = 0.1720. 


The numerical difference between the f2’s is 0.0190, and the numerical difference 
between the ji4’s is 0.0150. 

This example shows that Sheppard’s corrections are not valid to the precision 
to which they are usually given if they are to be used for the purpose of correcting 
raw moments. The last term in the fourth moment correction, 7/240, might 
equally well be, for example, —43/192 as in my special case. This will become 
more evident to the reader if he will draw the curve indicated in this example. 
To the eye it will appear exactly like the kind specified in the Euler-MacLaurin 
theorem; for example, much like the normal curve. Now suppose one adopted 
for the moment the point of view (which I have criticized earlier) of starting 
with the curve used in this example, breaking it up inte three partial areas and 
then finding the relation between the true and the raw moments. The partial 
areas found would be exactly those used in this example and this method would 
give us Sheppard’s corrections, but they would not be exactly correct, for in this 
instance my formulae give exactly the relationship between the true and the raw 
moments. The difference is due to the fact that in this instance the assump- 
tions permitting the use of the Euler-MacLaurin theorem in abbreviated form 
are not justified for this curve. But there is no way of telling at the outset, if 
one has given initially only the partial areas, whether precisely this curve or 
another which to the eye would appear very much like it is truly the curve which 
will graduate the same material when subjected to a finer classification. 





THE POINT BINOMIAL AND PROBABILITY PAPER 
By Frank H. Byron! 


1. An approximation to the sum of a number of consecutive terms of the point 
binomial may be found graphically and quite expeditiously by means of so- 


called “probability paper.’’ This paper is ruled so that the (2, y) graph of the 
equation of the integral of the normal curve 


1 * 5" 
y = —— é 2 dx (1) 
V/ 20 J—« 
is a straight line. Let the successive terms of the point binomial be represented 
as follows: 


(p+ 9)” =U tutes $Ueptrss +h, (2) 
where u; = .C.p”‘g' and p =q. Then the (z, y) graph of the equation, 


i.e., of the sum of first (¢ + 1) terms of this point binomial, is, in all but extreme 
cases, a set of points lying on a gently turning curve, so gently that its form may 
be represented closely by two straight lines, each passing through the median 
point as will be explained in the next section. As paper of this sort is readily 
obtainable, and as this method yields as great accuracy as is really useful in 
many problems, it is suggested that its use ought to be quite general. 


2. Sheppard’s Corrections. The formulae for the moments of the point 
binomial, mean = qn, o? = pqn, are exact without any corrections such as are 
used for grouped material. This fact has led us all (apparently) to assume that 
in fitting the curve to the point binomial one would get a better fit by equating 
the moments of the curve to the uncorrected moments of the point binomial 
rather than to the corrected moments. The studies made in connection with 
the preparation of this paper show that when the purpose is to equate areas to 
sums of terms the corrected moments should be used. The theoretical basis 
for this conclusion is as follows: 


To simplify the argument let us suppose that one were seeking that curve of 
Charlier type, 


F(x) = cobo(z) + cidi(z) + +++ Caga(x), (4) 


1 With the assistance of Burton H. Camp. 
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(where ¢, is the normal curve and q;, ¢2, --- its suecessive derivatives) whose 
integral would best fit the graph of (3). Since fitting is required only at the 
isolated points x = 3, 15, 23, --- , it is clear that one might obtain this by the 
two following steps. First let f(x) be any function whose integral meets exactly 
the requirement at these isolated points. What values this integral has at other 
points does not for the moment concern us. There are an infinite number of . 
such f(z) curves. Next let the c’s of (4) be so chosen that F(x) will fit f(x) as 
nearly as possible. The ordinary derivation of the c’s supposes that the fit 
between f(z) and F(x) is to be made by least squares, the residuals being weighted 
by the factor 1/+/¢(x). No matter what f(x) is chosen, the c’s can be deter- 
mined so that the weighted integral of (f(x) — F(x))? will be a minimum, but the 
value of this minimum will vary from one f(x) to another. We now desire to 
select that f(z) which will make this minimum value as small as possible, and 
it is reasonable to suppose that our best selection will be some f(x) which is as 
kindred to the nature of F(x) as possible. We shall not therefore choose an 
f(x) which oscillates wildly between the points where perfect fitting is required, 
(Fig. 1) nor yet an f(x) which is made up of the top bases of the point binomial 


afl ef 


Fic. 1 Fic. 2 


histogram; we shall prefer a modification (Fig. 2) of that histogram by a smooth- 
ing process. Such an f(z) will not have the exact moments of the point binomial, 
but, more nearly, those moments corrected for grouping. Then the determina- 
tion of the c’s will come out in terms of these corrected moments, not in terms of 
the uncorrected moments. (In fact the uncorrected moments would be the 
exact moments of an f(z) having an oscillatory character between the important 
points.) 

Of course, when 7 is large, the,difference is too small to be noticed and the use 
of Sheppard’s corrections is not worth while, and since n usually is large when 
approximations of this sort are needed, the point is not usually important. It 
was important in the computation of the tables of §4. Moreover, the use of 
Sheppard’s corrections does not invariably yield better results, the gain being 
masked sometimes by other effects to be considered in §3. An excellent illus- 
tration of uniformly better results is in fitting (} + 3)* by a curve of Type 4. 
The errors in the sums as derived from (4) with and without the corrections, is 
given on the following page. 
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t 2 3 6 7 8 


With 0001|—.0003!—.0001| .0000| .0001! .0903!— .0001| — .o002! 
Corrections | | | 
| ~ — | |- — 
Without | a on ~ 0000) — 0039) — .0022! — .0007| — .0001 
Corrections ie | Pe | | | 


3. The Stubby End. The other effects which mask this improvement are 
especially noticeable at the stubby end of a point binomial. We have to keep 
in mind here that the approximating curve (such as (4)), is required to turn a 
sharp corner, for, due to the least square method of fitting, it is just as important 
that it be close to zero when ¢ is negative, as it is that it be close to up, u,, --- 
when tis positive. Therefore, in order to turn this corner it has to dip below the 
x-axis in the neighborhood of f = — 3. Thismakes the approximating curve too 
low just to the right of ¢ = —4, unless the whole curve be arbitrarily widened. 
This arbitrary widening is customarily performed by not using Sheppard’s 
correction for o, and the result is a betterment of the fit at these points but a 
corresponding loss over the rest of the infinite interval. A good example? is 
(3 + 3)*. The fit is worse at the left end when Sheppard’s corrections are used 
but better over the rest of the interval. 

The same difficulty arises in another connection. If we compare the closeness 
of fit to a point binomial made by F(z) as written in (4) and by F(z) as it would 
be written if cy were zero, it often happens (as is well known) that the latter is 
actually slightly better on the average. How can this be true if the c’s are 
chosen by the method of least squares and the best choice as thus indicated 
makes c, different from zero? The answer is that the c’s are chosen so that the 
fit is best over the infinite interval, not merely over the interval from t = —} 
tot = n+ 3, and that furthermore the distant points are weighted more heavily 
than those near the center. Thus it might happen that a choice, other than the 
least square choice, and one in which cy would be zero, might be better for the 
restricted interval covered by the point binomial. This does happen especially 
when due to the abruptness of the stubby end of a very skew binomial, the 
curve has to dip below the axis in order to get by a sharp corner. A good ex- 
ample is the problem considered by Fry: (74 + 75). All the effects men- 
tioned are present here. The fit is on the average a little worse if c, is not equal 
to zero over the point binomial interval, a little better over the infinite interval. 


4. For graphical purposes a sufficiently good approximation to the median of 
(p + q)", is given by 


M = nq — (p — q)/6. 
2 The true values are given on page 220 of Mathematical Part of Elementary Statistics, 


by Camp, D. C. Heath and Company, 1931. 
3T. C. Fry, Probability and its Engineering Uses, p. 258, Van Nostrand, 1928. 
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The following tables enable us to find the first quartile Q;, and the ninth decile 
Ds. The accuracy to which they can be plotted is only about one-tenth that to 
which they are given here. Therefore accurate interpolation is seldom neces- 
sary. The values of S;,; are to be read from the graph at the points ¢ + 3, as 
indicated in the directions preceding the tables. The graphical method will be 
found efficient if one uses common sense in the computation. Numbers which 
are to be plotted should not be computed to a higher degree of accuracy than 
can be used graphically. In reading the values of S;4; it is well to remember 
that the true values lie on a curve, and that outside the interval from Q; to Ds, 
they are slightly less than those given by the straight line. Once the graph has 


9 
8 
7 
6 
5 
4 
3 
2 


345 6 7 6 4 95 
Yr axts 


Fig. 3 


been made, all the values of S,.; can be read quickly; it is not necessary to make 
a separate computation for each t. This method is therefore specially advan- 
tageous when one wishes to find several sums of this sort for the same point 
binomial. It should also be noticed that one can tell from the appearance of 
the graph about how far the true sum would be from the two straight lines and 
so estimate the error to which his reading is liable. 


5. Illustration. Find the sum of the first 7 terms of (3 + 4)”. 
Here ¢ = 6, M = 8.278, Q: = 6.726, Ds = 11.369. The graph shows that 
t 


>> = 0.224. The true value is 0.222. So the error is 0.002. 
0 
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An idea of the accuracy of the method is given by the errors (out of two places) 
that would be obtained for this point binomial for various values of t, as follows: 


¢ | 2] 4] 6 | 8 | 0 | 2 | 4 | 16 





Errors | .00 | .O1 | .00 | 00 | .00 | .00 | .00 | .00 


DIRECTIONS FOR USE OF THE TABLEs: Let p = q, M = ng — (p — q)/6, 
Q: = 1 + qn, Do = x2 + qn. On the graph draw the lines MQ; and MD,. 


Read S,4; at ¢ + 3. 


Values of x 


500 400 300 
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Values of x2 


400 300 75 





1.344 1.356 1. 1.481 See Auxiliary 
323 330 396 Tables 
314 319 367 


309 313 352 
305 309 342 
303 306 335 


301 304 329 
299 302 "325 
298 300 321 


296 299 318 
297 313 
294 308 


290 291 301 
288 289 297 
286 287 288 293 


285 286 286 290 
284 284 285 287 


282 282 282 
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INEQUALITIES AMONG AVERAGES 


By Nirvan Norris 


Numerous inequalities among averages of various types are condensed in the 
monotonic character of the function 


1 


¢(t) = (2 +2%34--- + #2) 


n 





of the positive numbers 2, %2, --- , Xn, not all equal each to each. Fort = —1 
this function is the harmonic mean; for ¢ = 0 it is the geometric mean; for ¢ = 1 
the arithmetic mean; and for t = 2 the root mean square. The relations 
among these four means which customarily are proved by special and dis- 
connected methods appear easily as applications of the theorem that ¢(t) is 
an increasing function of t. Thatis, for any values of t; and & such that — 
<t<#< + ~, it will be true that (4) < (tf). Several proofs of this theo- 
rem have been published, many of them very complex. An extremely simple 
proof is herewith presented. 

That ¢(t), ¢’(t) and ¢’’(t) all exist and are continuous for all real values of t 
may be shown by expanding each of the quantities x; in a series of powers of ¢ 
and considering the remainders after each of the first three terms. The ordinary 
rule for evaluating forms reducing to 0/0, which requires the function under 
consideration to be continuous and to have at least a continuous first derivative 
for t = 0, may then be applied to [log ¢(¢)]/t to show that $(0) is the geometric 
mean. It is clear that ¢(— «<) and ¢(+ ~) are respectively the least and the 
greatest of the z;. This fact and the monotonic property of ¢(¢) make it evident 
that for each real value of ¢, the function may be regarded as an average in the 
usual sense that it lies within the range of the observations. 

For a simple demonstration of the increasing character of $(¢), consider the 
auxiliary function 


— P(t) _ si ze} Setlogr | Eat 
F(t) =t $(t) =t a log t —_—-— — log ; 


t wt nr 


It is clear that ¢’(t) has the same sign as F(t). The theorem will be proved by 
showing that the sign of F(t) is positive for all values of ¢ except zero, when 
¢’(t) vanishes. 


1 Professor Harold Hotelling rendered invaluable assistance in condensing for publica- 
tion the material herein presented from a more extended study of generalized mean value 
functions. 
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Differentiating the last expression with respect to ¢t, one obtains upon sim- 
plification 





aan 


- F(t) = 2 


[(S2') (2 xt log? xz) — (Zz log z)?]. 

By Cauchy’s inequality (known as Schwarz’ inequality when applied to integrals 
instead of sums), the expression in square brackets is positive. Hence F’(t) 
has the same sign as ¢t. Consequently F(t), since it diminishes for negative 
values of ¢ and increases for positive values, has a minimum fort = 0. But by 
direct substitution, F(0) = 0. It follows that F(t) and ¢’(t) are positive for all 
values of ¢other than zero. Therefore ¢(¢) is an increasing function. 

By direct general methods it is possible to show that 


¢'(0) = (Iz) = [n>(log x)? — (3 log x)?]. 


2 
This expression obviously vanishes only when nz (log x)? = (= log x)?, a condition 
which is satisfied only in the trivial case when x; = x2 = --- = Xp. 
A proof exactly parallel to that given above may be applied to integrals or, 
more generally, to Stieltjes integrals. The monotonic increasing character of 


=0 
function integrable in the Riemann-Stieltjes sense, such that y(o) — ¥(0) = 1, 


© 


and such that i x‘d(x) exists for every real value of t. In terms of statistical 


z=0 


* 1 
| [ wdy(z) | appears in this way if one assumes that ¥(x) is a non-decreasing 


theory, this consideration extends the theorem from samples to populations of a 
very general character. 

Proof of the increasing character of ¢(t) has also been derived from Hdélder’s 
inequality, the demonstration being expressed in terms of Stieltjes integrals.? 
The simplest general proof of the monotonic attribute of ¢(¢) heretofore published 


appears to be that of Paul Lévy.* As early as 1840 Bienaymé‘ presented a 
generalized form of ¢(t), namely, 


1 
(18 + cost 2 aiid tet 
Ct+eeters + en 


and announced, without proof, its increasing character. In 1858 a proof of the 
monotonic quality of ¢(t) for special cases was published by Schlémilch.’ Of 


? 





2 J. Shohat, “Stieltjes Integrals in Mathematical Statistics,’’ Annals of Mathematical 
Statistics (American Statistical Association, Ann Arbor, 1930), Vol. 1, No. 1, p. 84. 

3 Calcul des Probabilités (Gauthier-Villars et Cie., Paris, 1925), pp. 157 f. 

4 Jules Bienaymé, Société Philomatique de Paris, Extraits des Procés-Verbaux des Seances 
Pedant L’Anée 1840 (Imprimerie D’A. René et Cie., Paris, 1841), Seance du 13 juin 1840, 
p. 68. 

5 Q. Schlémilch, ‘‘Ueber Mittelgréssen verschiedener Ordnungen,”’ Zeitschrift fiir Mathe- 
matik und Physik (B. G. Teubner, Leipzig, 1858), Vol. 3, pp. 303 f. 
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the more recent general proofs of the increasing character of ¢(¢) which have 
appeared, those of Jensen,’ Pélya,’ Jessen,’ and Carathéodory® may be men- 


tioned. A recent application of ¢(¢) to index number theory is that of Professor 
John B. Canning.” 


VassaR COLLEGE. 


6 J. L. W. V. Jensen, ‘‘Sur Les Fonctions Convexes Et Les Inegalités Entre Les Valeurs 
Moyennes,’’ Acta Mathematica (Beijers Bokférlagsaktielbolag, Stockholm, 1905), Vol. 30, 
pp. 183-185. 

7G. Pélya and G. Szegé, Aufgaben und Lehrsdtze Aus Der Analysis (Julius Springer, 
Berlin, 1925), Vol. I, pp. 54 f. and 210. 

8 Borge Jessen, ‘‘Bemaerkninger om koveskse Funktioner og Uligheder imellem Middel- 
vaerdier,’’ Matematisk Tidsskrift (Charles Johansens Bogtrykkeri, Copenhagen, 1931), 
No. 2, 1931, pp. 26-28. 

9 Attributed to Professor Constantin Carathécdory in an unpublished manuscript of 
Professor Harold Hotelling. 

10 “A Theorem Concerning a Certain Family of Averages of a Certain Type of Frequency 
Distribution,’ a paper presented before a joint meeting of the American Statistical Asso- 
ciation and the Econometric Society at Berkeley, California, June 22, 1934. 





































MATHEMATICAL EXPECTATION OF PRODUCT MOMENTS OF SAM- 
PLES DRAWN FROM A SET OF INFINITE POPULATIONS 


By Hyman M. FEtLpmMan! 


Introduction 


In the second part of his investigations, “On the Mathematical Expectation of 
Moments of Frequency Distributions,” Tchouproff presented a method which 
may be interpreted as sampling from a set of infinite univariate populations. 
In the present paper this method is extended to the study of moments of product 
moments of samples drawn from a set of infinite bivariate populations. It is 
also shown how this method may be extended to populations of higher order by 
deriving some of the simpler formulae for populations of three and four variables. 

Tchouproff’s method has been criticised* because of the complicated algebra. 
On close examination it is found, however, that it is not the algebra which is 
complicated but rather the symbolism. Tchouproff introduced a great variety 
of symbols both in his derivations and in his results. As a consequence his work 
seems very intricate. If, however, the number of symbols is reduced, and the 
symbols themselves are simplified, which can be easily accomplished, the under- 
lying idea of Tchouproff’s method is found to be very simple. 

Quite a complete study of product moments of any bivariate population has 
been made by Joseph Pepper in his “Studies in the Theory of Sampling.’* His 
method is essentially an extension of Church’s®’ method, in his studies of univa- 
riate populations, to bivariate populations. He does not, however, derive any 
generalized formulae. In the present study generalized formulae for both the 
first moment and the variance of product moments of any order are obtained. 

It may be noted here, that all of Pepper’s formulae for any infinite population 
can be obtained from those of the present study as special cases, by assuming 
that all the populations in the set are identical. 


1 A dissertation presented to the Board of Graduate Studies of Washington University in 
partial fulfilment of the requirements for the degree of Doctor of Philosophy, June 1933. 

2 Biometrika, Vol. X XI, Dec. 1929, pp. 231-258. 

3 Church, A. E.R. ‘On the Means and Squared Standard Deviations of Small Samples 
from any Population,’’ Biometrika, Vol. XVIII, Nov., 1926, pp. 321-394. 

4 Biometrika, Vol. XXI, Dec. 1929, pp. 231-258. 

5 Church, A. E. R., ‘‘On the Means and Squared Standard Deviations of Small Samples 
rom any Population,’”’ Biometrika, Vol. XVIII, Nov., 1926, pp. 321-394. 
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CuHapTerR I. Notations and Definitions 
Let (Xi, Y1), (Xe, Y2), --- (Xn, Yn) be n bivariate populations each following 


any law of distribution whatever. The product moment of order a in X and b 
in Y of the kt population will be denoted by P£,. It is defined as 


Pos = E(Xz — ax)* (Y; — by)” (1.11) 
where a, = E(Xx), b, = E(Y;,), (1.12) 


and where the symbol E£ signifies the expected value or the mathematical expec- 
tation of a quantity. 

Regarding each of the » populations of the set as infinite, samples of n are 
drawn, each member of a sample from one of the n populations.’ The individual 
which is drawn from the k** population will be denoted by (2x;:, yx); and the 
product moment of order a in x and b in y, of such a sample will be denoted by 
Par. This product moment may then be defined as 


Pa = nN" S (x. — x)* (yx — y)? (1.13) 
where z= n Szx,, y = n Syx. (1.14) 
The symbols a and b will now be defined by the equations 


a= n Sax, b= n-} Sb;.. (1.15) 
Obviously E(x) = E(n Sa;x) = n SE(X;) = n Sa, = a. (1.16) 


Similarly E(y) = b. That is, the mathematical expectation of the mean, of 
such a sample as was described above, is equal to the average of the means of all 
the populations.*® 


In order to make the equations as compact as possible the following additional 
symbols will be employed: 


Le — Ax = Uk, r—a=u, and uz, — u = U;, 
(1.17) 
yx — by = UE, y—b=2, and vy. — v = V; 
also a;.—-a= Ax, b;. _— b = B,.. 
From the above definitions it easily follows that 


E(ux) = E(vx) = E(U,) = E(V;) = E(u) = Ev) = 0. ~— (1.18) 


6 The term infinite is used here in the probability sense. It is defined very clearly by 
Church in his ‘‘Means and Squared Standard Deviations of Small Samples,’’ Biometrika, 
Vol. XVIII, Nov., 1926, p. 322. 

7 It may be easily shown that this is equivalent to drawing a sample of n from a set of 
any finite number of populations. The number drawn from each population, however, 
must be specified. See Biometrika, Vol. XIII, 1920-21, p. 295, footnote. 

8 This, of course, is a result of the Lexis Theory, for Poisson and Lexis Series. 
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The notation is now completed with the definition of the symbol Q;; by the 
equation: 


























Q;; = S(a, — a)é(b, — b)i = SAiB). (1.19) 








CuapTer II. The Mathematical Expectation of p., 





The mathematical expectation of pa» will be denoted by j.s. In the terminol- 


ogy of moments this would be called the mean or first moment of the distribution 
of Pad. 





1. The Mathematical Expectation of pi. 
the expected value of puis ju. By definition 


According to the above notation 





Pu = E(pu) = EnS(z; — x)(yi — y), (2.11) 





and obviously En“S(a; — x)(yi - y) = n“SE(xi — x) (yi — y). 
Writing 


a; — x = [(x; — ai) — (e — a)] + [a — a] 
ys —y = [yi — bi) — YY — BD] + [Bi — 4] 
equation (2.11) may be written as 
pu = n-SE(Us + AQ(Vs + Bd 
= nSE(U.V;) + n“SAiE(V;i) + n“SB;E(U;) + nSE(A;B)). 


U;+ A; 
V; + B;, 


I 


Since for any given population A; and B; are constants, it follows that 
E(A iB;) =A iB. Hence 


n=SE(A iB i) = nSA iB; = n On ‘ 





Making use of (1.18), it is seen that the terms n“SA ;E(V;) and n“'SB,E(U,) 
are zero. The only term left to evaluate is therefore n“'SE(U;V;). Since U; 
and V; are symmetric functions of the corresponding small letters, their product 
is symmetric in uw; There is therefore no loss in generality if attention is 
concentrated on a single subscript, say 1. 

We may therefore write 


n 'SE(U;Vi) = n“E(U,V;) + nSE(U;V)). 
“2 
Remembering that U; = u; — u = ui; — nSu;, we may write, 
U; =u; —-—u =u; — nu + uw + --- + u,) 
=n [nui — (wm +m tee) $+ Wat Winn +--+ + u,)] 


*The 2 at the bottom of the S simply indicates that the summation begins with 7 = 2, 
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where m; = n — 1. In general, n; will denote the number n — 7. Similarly 
Vi = nmin — (01 + 02 + +++ + Osa + Vi4n + +++ + 2). 
Thus 
nSE(U Vi) = nF (mus — Ue — +++ — Un) (MW — Ve — +++ — Vp) 


+ n3SE(nu; — Um eee — Ue — U4 — <-> Un) 
2 


(nv; — Vy — +++ — Vi — Vig — +++ Un). 


When the right hand side of the last equation is expanded the only terms which 
appear are of the form E(u,;) and E(uw;). The last one must vanish for u; 
and v; are independent and hence E(ujw;) = E(u;)E(v;) = 0. From the last 
equation above it is easily seen that the coefficient of E(u) is 


n(n? + mn) = n3n(m + 1) = n?n;; 


and because of the symmetry this is obviously the coefficient of any term of 
that form. Hence 


nSE(U Vi) = nn, SE(ui) . 
Since u; = 2%; — ai, vi = yi — 5;, then 
E(uw,) = E(x; — ai)(yi — bs) = E(X; — a;)(¥; — bs) = Pi, 
and in general, 
E(ujvj) = Pi ;. 
We thus get the formula 
pu = n—n,SP}, + nQn. (1) 


Now suppose all the n populations are identical. Then all the A’s and also 
all the B’s vanish and therefore, Qi, = 0. The formula (1) thus becomes 


| ‘ 
pu = Py. (1’) 


This is exactly Pepper’s formula for pu for an infinite population.® 


2. The Mathematical Expectation of p.. By definition 


pu = EnS(x; — x)*(yi — y). (2.21) 


® Biometrika, Vol. XXI, p. 233, Eq. A, N = ©. As was already stated in the introduc- 
tion, all the formulae of the present study reduce to Pepper’s when the above assumption 
is made. 
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Proceeding as above it is seen that 
En“S(x; — x)*(yi — y) = n“SE(xi — 2)*(yi — y) 
= nSE(U; + Ai)*(Vi + Bi) = n“SE(UZV,) + 2n—SE(U;ViA,) 

+ n“SE(U?B;) + n—SE(V;A2) + 2n1SE(A,BU,) + n“SE(A{Bj) - + + - (2.22) 


It is quite evident that the two terms before the last vanish. To evaluate the 
remaining terms, we employ the reasoning of section 1 of this chapter and write: 


SE(U2V,) = E(U2V;) + SE(U2V)) 
2 





= nE (niu — Ue — --- )(md — v2 — ---) + nSE(mu; — um — ---) 
2 


(nvi — 4 — <-->), 


Since terms of the form E(u?v;) vanish, only the coefficient of the term E(u7v,) 
must be found. Again considering the subscript 1, the coefficient of E(uj1) is 
easily found from the last equation to be 
n(n? — m) = n-n(my + 1)(m — 1) = n-ne. 


Thus 









nSE(U2 V;) = n-nyn,SE(u? v;) = n—nyn2SP3, . 





For the second term of (2.22) we have 





SE(UV;:A) = E(UiV;A;) + SE(U.V:A)) 


= nF (nyu — Ue — --- )(mvy — ve — --+ )Ar + n?SE(mui — w — ---) 


)A é. 





(uv; — % — <--> 





The coefficient of E(u,v;) in the first term of the right hand side of the last 


equation is n~*n{A;. In the second term it is n~*SA ; = —n~?A,, since SA; = 0. 
2 








It therefore follows that 





2nSE(U ;V;A i) = 2n-*noSPi 1A i. (2.24) 


Quite similarly 





n—SE(U2B;) = n-2n2SP3 Bi, (2.25) 






and it is obvious that 





n—SE(A2B;) = nQz. (2.26) 






* Note that the u which has the coefficient n; does not occur among the u’s which have 
the negative sign. 
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We thus get the formula 
pa = nn n2SP 3, + n-ne S(2Pj, A; + P3B;) + nQn. (2) 


3. The Mathematical Expectation of ps, and po. 
pu = EnS(a; — 2)*(yi — y) = n“*SE(x; — 2)*(ys — y) 
= nSE(U; * A;)*(V; + B,) = nS{E(UiV, + USB; + 3U7V,A, 
+ 3U7A,;B; + 3U;V;A? + 3U,A7B;.4+ V;A’* + AB}. (2.31) 
The two terms before the last are zero. The last term is 
n—SE(A?B;) = n—Qy. (2.32) 
By (2.23) and (2.24) and some slight manipulation 
3n-1 SE(U7A;B; + U;V;A7) 
= 3n-n,S(P},A;B; + Pi, A?) + 3n(Q,,SPi, + Qi,SPi,), -(2.38) 
and by (2.22) 
nSE(UiB,; + 3U2?V,;A,) = n(n? + 1)S(Pi,B; + 3Pi,A). (2.34) 


The only new term which is to be evaluated is SE(U?V;). This may be 
written as follows: 


SE(U?V,) = n“*SE(mui — um — --- (mo; — 1 — ---). 


When the right hand side is expanded it is found that the only non-vanishing 
terms are of the form E(U?V,) and E(u?u,v;). Only two subscripts, therefore, 
have to be considered. Without any loss in generality these may be taken as 
1 and 2, and the right hand side of the last equation may then be written as 
follows: 


SE(nyu; — uw — --+ (ny; — 4 — +--+ ) = Eins — ue — --- )F(nyv — 02 — -- 


+ E(nyu2 — us — ---)8(my1 — ve — ---) + SE (mui — us — Ue —---)3 
3 


° (nv; — 11 — V2 — +++), 


From this last expansion it is easily seen that the coefficient of E(uiv,) is (n{ + n,) 
and that of E(u}ujv;), (6n? + 3n) = 3(2ni + m). We thus finally obtain 


SE(U;V;) = n-{(nj + n,)SE(uj v;) + 3(2n} + )SE(uju;;)}. 


But by (2.12) E(ujv,) = P3,, and since u, and u; and u; and v; are independent 
E(uzu;v;) = E(uz)E(uj;v;) = P3>Pi,. Whence 


E(U3V,) = n-{(nt + n,)SPi, + 3(2n? + m)SPi,Pi,}. (2.35) 
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From (2.31) and the succeeding equations we finally get 
pu = n-*{(ny + m)SP3, + 3(2nt + m)SP*20P ii} 
+ n“{(n} + IS(P3 0B; + 3P2,4,)} + 3n{ (nj — 1)S(P2 04; B; + Pi1 47) 
+ Qu SP 20 + QoSPii} + 27Qa . (3) 


The derivation of pz is so similar to that of jz, that it would be mere repetition 
to go through the details again. We shall therefore merely write down the 
formula for p22 which is 


p2 = n—{(n} + n,)SP3 > + (2ni + n)S(P 30 be + 4P;,Pi,)} 
+ 2n-{ (nj + 1)S(P3 ,B; + P;2A;)} + n-3{ (nj = 1)S(P;B; + 4P;,A;B; 
+ P5245) + QoSPo2 + Qe SP20 + 4QnSPii} + n*Qee. (4) 





4. The Mathematical Expectation of the General Product Moment p,». 
So far, formulae for the mathematical expectation of ps, for particular values 
of a and b, have been derived. The method used in deriving these is, however, 
perfectly general, and now, that it has been sufficiently illustrated, it can be 
easily generalized. 

By definition we have 


Des = Eln-'S(x; — 2)*(ys — y). 


Making use of the notation of Chapter I this may be written as 
a,b n 
npo = ES(U; + Ai)(Vi+ Bi)’ = S CeC?SE(US*V}-"A%B?) (2.41) 
q,7r=0 1 


where 
a! b! 
~ gi(a —q)!’ ri(b—r)!" 
Expressing the U’s and V’s in terms of the w’s and v’s and setting a — q = l, 
b — r = m; we may write for a particular pair of values qg and r: 


nitm SE(U! V"™A‘B‘*) = SE(nu; — um — ---)'(nw; —, — ---)™AYBi. (2.42) 


c C* = 


Consider, now, the general term in the expansion of the right hand side of 
(2.42). Itis of the form: 







I'm! 


Ila, !11B;! (— 1)'*™(—n,) 77 E(nu5} Sahte uskyer eee vik A1B*) ; (2.43) 


74 









where ITa,! = a;!ae! +--+ a,x! 


*In this case, and also in the formulae that follow, whenever two or more indices 
appear in a summation, it will be understood that no two of them can have the same 
vajue simultaneously. 
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For particular sets of values ji, jo, --- jx, a1, @2, +++ ax, and Bi, Bo, --- Bx, this 
term will appear in every member of the summation of the right hand side of 
(2.42), and its coefficient will differ only in the exponent of (—7,) and in the 
subscript 7 of A?B’. Because of the symmetry there is no loss in generality if 
we take for ji, jo, --- jx, the first k integers. We now break up the summation 
of the right hand side of (2.42) as follows: 


SE(mu; — um — ---)'(mv; — v, — ---)"ALB; 
1 


= E(mum — uw — ---)'(mv, — vo — ---)"ATBi 
+ E(t, — uw — ---)'(niv2 — 0) — ---)"AGB3 +--+ + (mu, — w — - 
(n% —%—---)"ATBi + S Elnu;— wm — ---)' 
i=k+1 
(mvj — v1 — ---)"ATB?. (2.44) 


From (2.44) we easily get for the total coefficient (excluding the numerical 
factor) the expression 


k n 
S (—n,)*"*A, Bi + S AjB;. 
h=1 h=k+1 
Writing 


&+14 


n n k k 
S Aj B;, = SA{B;, — SA{B;, = Q,, — SALB;, 
1 1 1 


the general term, (2.43), together with the total coefficient, may then be written 
as 


ltm! f* h 
(—1)tm £2 ™ =) og (norton — 1] ATBL 4+. Q,,5 E uttot, 
‘ h=1 


Ila, ! IT), ! 


h= 


Since u; and u;, v; and v;, and u; and v; are independent while u; and v; are 
not, we have: 

I. Ellu, »» = Eu, v, = UP, 5, 

II. Any term in which a, + 8, = 1 must vanish. 

From II it follows that the maximum number of subscripts which can appear 
in any term in the expansion of (2.42), i.e. the upper limit of k, which will be 
denoted by t, cannot exceed (J+ m)/2. In fact when! + mis even, t = (1+ m)/2, 
while when | + m is odd, t is the largest integer less than (1 + m)/2. 

Making use of (2.41), the equations following it, and the reasoning of the last 
paragraph, we finally get the formula: 

r e-g.b—-r ft 
(=n) 8" 


n a,b 
n(—n)*+*p. = (a!) (b!) S S S 


ja=1 q,r=0 q! r! an=0,8,=0 k=1 


: ih 
{ S [(—m)+6 — 1] A$, Bi, + Qu 288A 
1 


h= an! Br! 
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The following restrictions on the a’s and #’s must be observed 


(a) ar +a+-+---+a,=a-—q 
(bt) i+ h+---+6.=b-—47 
(c) on +P #1. 


In case the n populations are identical (5) reduces as follows: For q 
r= 0, Af = 1, Bf =1, and Qu = n; while in every other case A%B, 
Q,, = 0. The summations with respect to qg and r, therefore disappear. 

Consider now the summations 


s Pi P22... Pik , 
s2"1. gi 1 
Since all the populations are the same we may drop the j by actually carrying 
out the indicated summations. If, then, there are c repetitions among the k 
pairs of integers apBn, in which a18;, arBe, --- a-B-, are repeated 1,, 12, --- 1, 
times respectively, then we have; 


k! Cy 
wee PS ee gee Pie - 
ae sank Bh Lth!---b! hBh 
We thus arrive at the following corollary: The mathematical expectation, 
Pav, of the product moment, pas, in samples of n from a single infinite population 
having any law of distribution is given by 


a,b ! ! i n af n 
n(—n)**pap = S (a ') (b ') Sg | S (— 1) 2rtBr + ns | on TI Pons (5’) 


ouprwe Ties | TBs Seo: -1,! 


h=1 


Note: In deriving these general formulae it was assumed that n > t. There 
is however, no loss in generality in this assumption. For, if t > n, we may 
suppose that, tn41 = In42 = --- = 2% = 0, and hence P23’ = --. = Pj, =0, 
and thus the above reasoning is still valid. 


5. Formulae for pi, ps2, ps1, P12, P33. Formulae for pa, in which a + b = 5, 6, 
7,8 have been obtained. But for (a + b) > 6 these formulae become very long, 
and since these will be of no use-in the subsequent work, only those of order 5 
and 6 are given below. 


+ n-5 {(nj a n)S(Pi B; + 4P:,A,) + 6nnS(P3, BP, + 2P''A;Pi,)} 
* This is a generalization of Pepper’s results for N = ~. See Biometrika Vol. XXI, 
pp. 231-240. 


+ The symbol P},4;P3o is an abbreviation of the full term (A; + A;) (Pi1P40 + 
P{,P30). Similar abbreviations will be used in the other formulae. 
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+ 2n- {(n} + 1) S(2P3,A:B; + 3P2,Aj) — 2QuSP3_. — 3Qn0SP3,} 

+ 2n-* {(nt — 1)S(2P},A$3P3,A7B,) + 2QsSPi, + 3QxSP3o} + nQu. (6) 
po =n {(n} — m)SP32 + nn3S(P3oP 2 + 6P31Pi, + 3P20Pis)} 

4+ n- {(ni — 1)S(2P3,B: + 3Pi.A;) + 3nnS(Pi, Pi, B: + (PioPi: 

+ 4P1,P{i) Ad}+ vn (my + YS(P3 Bi + 6P2, AB; 

+ 3P',A2) — QuSPi, — 6QuSPi, — 3Q0SPi,} + n-* {(n? — 1)S(3Pi,A;B? 
+ 6Pu A.B; + Pp 2A¢) + 3QuSP 35 + 6QxSP}, + QsSP52}+ n-"Qze . (7) 
pa =n {(n§ + m)SP§, + 5(ni +n? + n)S(PioPi, 

4+ 2P3, Pi.) — 10(2nj — 2) SP}, Pi, + 3Gn? + ndSP;,Pi.Pi.} + wo *{at 
+ I)S(P5 0B: + 5P41As) + 10(my + 1)S*2P3 oP 2 0B: + (2P30P iy 

+ 3P31P}o)Ail — 10nneS*[2P5 Pi 9B; + (2P3oPi1 + 3P21P2o)Aj)} 

+ 5n- {(n{ — 1)S(Pi,A:B; + 2P31 47) + 6nnS(P3 oP} AB: + 2P3oPi,47) 
+ QuS(Pio + 6P30P30) + 2Qx0S(P3, + 6P3oPi,)} + 10n-*{(n} 

+ 1)S(P3,A{B;+ P3},A%) — QaSP§ 5 — QSP3,} + 5n-{(n? — 1)S(2P3,A3B; 
+ P{,Aj) + 2QsSP 29 + QwSPi1} + 2-Qn . (8) 
po = n{(ny + m)SPis + (nt +i +m)S(PioPie + 8P3iP ii 

+ 6P32P3o) + 4(2n} — m2) S(P3oPi2 + 3P21P21) + 63nj + n3)S(P2oP2 oP o> 
+ 4P3oPiiPi1)} + 2n{(ni + 1)S(Pi 1B; + 2P32A,) + Anz + 1)S[(2P3 oP} 
+ 3P§ Pi.) B: + (PiePic + OP3,P/{, + 3Pi,PidAd — neS(2P;, Pi, 

+ 8P2 P30) By +(P3oPi2 +6P21Pi1 +3Pi2Pio)Ai} +n {(mi — YS(Pi Bi 
+ 8Pi, A,B; + 6Pi,A?) + 6nn,S[Pi, Pi ,B? + 4P}, Pi, AB: + (Pi,Pi; 

+ 4P},P{,)Ai) + QuS(Pio + 6P20P20) + 8QuS(P3, + 3P30P}1) 

+ 6PoS(P22 + P2oPi2 + 4PiiPii)} + 4n“{ (ni + 1I)S(P3 AB; 

+ 38P3,A7B: + Pji2Aj) — QuoSP3o — 3QxnSP2, + Q30SPi2} 

+ n{S(6P3,A7Bi + 8P{,AiB:i + Po2Aj) + QuSPo2 + 8QnSP}, 

+ 6Q22SP 3} + Qe . (9) 


Bos = n-"{(n§ + m)SP33 + 3(ni + nj + m)S(P3,Pi2 + 3P32Pi, 
+ Pi 3Pi.) = (2n} — N) S(P3o is + 9P3, 9) + 9(3n} + n3)S(P3oPi, P52 


*The repetition of this expression signifies that A and B factors are coupled only with 
those P factors which have corresponding indices. 
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+ 4Pi,Pi,Pi1)} + 304 (mt + YS(P32Bi + P2sAi) + (ni + 1D) Sl(P30P ie 
+ 6P2, Pi, + 3Pi2P20)B:i + (PosP20 + 6Pi2Pi, + 3P21Pi2)Ail 

— nneS[(P3oP52 + 6P2,Pi, + 3Pi2P3o)B; + (PosP2o + 6Pi2Pis 

+ 3P2,P)2)Aj]} + 3n{(nj — 1)S(P3, Bi + 3P2,A:B; + P} 343) 

+ 3nyn2S[P2 P11 Bi + (P2oPi2 + 4Pi,P}1)ABi + Pi1Pi2A47)) 

+ S[Qee(P31 + 3P30Pi1) + 3Qu(P22 + P2oPi2 + 4Pii:Pi1) + Qu(Pis 

+ 3Po2Pi1)]} + n*{(n} + DS(P3,.Bi + 9P2,ABi + 9P{.A7A: + Po3A}) 
— S(QusP3o + 9Q:2P31 + 9QxP is + QP53)} + 3n-{ (ni — 1)S(P2 AB 

+ 3P;, Aj Bi + Po,A{Bi) + S(QisP20 + 3Qx2P i, + QaiPi2)} + Qss. (10) 
CuHapTeR III. The Mathematical Expectation of the Variance of p., 


1. The Symbols .mp,, and »M>,,. Denoting the variance of pa» by m and 
the mathematical expectation of 2mp,, by 2M», , we have the definition, 


2Mrg, = {n— S(x; — x)*(ys — y)® — Pav}? 





= n-?S%(ax; — x)*(y; — y)® — 2n paw S(z; — x)*(y; — y)® + p2,, and 
2M, = E(x) = E{n*S*(x; — x)*(yi — y) — 2n-PavS(x;i — x)°*(yi — y)® + Dav} 
= nE[S(x; — x)**(yi — y)™] + 2nE[S(x,; — x)*(2; — x)*(yi — y)*(y; — y)?I 
— 2n-parE[S(x; — x)*(yi — y)*] + Diy = NPrar 
+ 2nE[S(x; — x)*(yi — y)*(x; — x)*(y; — y)*] — piv. (3.11) 






Before attempting to expand the right hand side of (3.11) for any values a, b 
we shall derive the formula for »M>,, to illustrate the procedure. 










2. The Mathematical Expectation of 2m»,,. By (3.11) we have 





Mr, = Noo + 2n7E[S(x; — x) (yi — y) (a; — x)(y; — y)) — Bis (3.21) 













The first term is given by (4) and the last by (1). The only new term is the 
middle one. To expand it let us write it in terms of U and V. We then have: 


n~SEl(zi — x)(yi — y)(ai — z)(yi — y)] = n SEU: + Ai)(Vi 
+ Bi)(U; + Aj)(V; + B;)] = n?{SE[U V UV; + (UiViU;B; + U;V;U.B;) 
+ (UiViV;A; + U;V;ViA,) + (UiV:A;B; + U;V;A<B;) 

+ (UiV;A iB; op U;ViA iB;) + U U;B;B; ot ViV;A iA; a 4 vanishing terms 
+ A,B.A;B)]}. 








(3.22) 
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The evaluation of the last term is very simple. For 
SE(A;B;A;B;) = S(A,B;A;B)), 
and from the elementary theory of symmetric functions we have: 


2(A,B,) — S(A2B? 
S(A;B,A,B,) = ce See . 


2 _ S(A2B? _ 
SE(A,B,A,B,) = a - Sis = Se, (3.23) 


To expand the first term and also the remaining ones, we return to the u, », 
notation defined in Chapter I. We then write 
SE(U VV <U;V;) = n4SE[(mus — uw — --- (my — 1 — «> ) 


(muy — uw — +++ )(mvj — 1 — --- dI. 
The only terms which can appear in the expansion of the right hand side of the 
last equation have the following form: 


E(ujvi),  E(ujv;), E(upur;) , 


ae 


i.e., exactly those which appear in the evaluation of jx2. Remembering the 
symmetry, there will be no loss in generality if we take for 7 and j the integers 
1 and 2. To find the coefficients of the three characterstic terms, the above 
summation may be broken up as follows: 


nSE(U,V;U;V;) = El (mim = ae = eA -)(nw1 a a -)(nyu2 = +4 -) 


(nye — 4, — ---)) + E{[nim — we — ---)(mimi — » — +--+) + (rie — w 


= -)(nive — a = = -)|S(nu; =a = =. -) (nw; ty eee -)} + SE[(niu; 


—U-—-- -) (ny; —U—:- -) (niu; —uUW—-:: -)(nyo; —Ui— °° -)] ° (3.24) 


Writing the three terms in a row and their coefficients from the three parts of 
(3.24) in columns below these terms, we get the following scheme: 


E(ujvi) E(u2v5 + uiv}) E(uy v1 U2 v2) 
n? ni (n* + 1)? 
n(n? + 1) —2nin2 2n3 
N2N3 NeNs 


5 7 amas 


a» 2 n(nt + ni — 3n, +3), 
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With the aid of the above equations we finally get: 










* 
SE(UV.U,V,) = ae ee =P api, ~ Sari ri, 





+ n(nn? — 3n,)SP}, is} 
Proceeding in the same way we find: 
SE(U,V,U;B; + U;V;U,B,) = n-*(2n> + n,) SP}, B; 
SE(U,V,V;A; + UjV;V;iA) = n—(2n? + n,)SP},A; 
SE(U,V;A;B; + U;V;A;B) = —nn.SPj, A,B; + (n? + n,)Q;,SP}, 
SE(U,V,;A;B; + U;V;A;B;) = 2nSP}, — QySP}, 
SE(U,U,;B,B; + ViV;AiA;) = nS(P3,B? + P5.A2) — 3S8(QuP do + QuP20)- 
Collecting terms and simplifying we finally get: 


2M py, = n-*{n}SP 3 + S(P 26 be + 2P}, Pi) — n’S(P},)"} 











+ 2n-*n,{S(P2,B; + Pj2A,)} + 0 {S(P2oBi + 2P},A;B; + Po2Aj)}. (11) 






Corollary 1. 
(11) becomes 


In case X; = Yj, i.e., when the set of populations are univariate, 





Mog = n-*{n?S[Pi, — (Pi,)*] + 4SPi, Pig} + 4n-n, SPi,A,; + 4n—SPi,A?. 
(11’) 








This is Tchouproff’s formula for the expected value of the variance of samples 
of n.?° 


Corollary 2. In case the n populations are identical (11) becomes 





2M ry, = n-*ny[n; Poo + Poo Po: — P31). (11’”) 


3. The Mathematical Expectation of »VW»,,. We now return to the general 
equation : 






Mra, = "Pog, — Day t2n-? S E(x, —x)*(y; — y)*(2; — z)*(y; — y). (3.11) 


¢*1,j371 


* Since E(uzv?) = Pe, E(u7v?) = P30Pis, etc. 
10 See Biometrika, Vol. XIII p. 295. 
1! See Biometrika, Vol. X XI p. 234, Cor. 1. 
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The first two terms are given by (5). To evaluate the last term we write: 


SE|(x; — x)*(ys — y)*(a; — x)*(y; — y)*] = SE|(U; + As)*(Vi + By)*(U; + A;)2 
a,a,b,b, 
(Vi + B)'] = SE(USViU5V5) + SOF Cr Cr Cr, 


ryT ar ger 4 =O 
SE(USV1U8V2A1BT3A 2B") = 20) SE{ (niu; — «+ +)2(nyw; — - ++) 
1 


. -)*(nyv; ie .)e oo S mlritretr3trs) C., ene C?, SE|(nyu; = 06 .)@ 


rs 


— ---)(nyu;j — ---)4(nyw; — -- -)®AIBi3A 2B, 


wherea =a—1,8B=a—ny=b—7r3,6=b— 7. 


The right hand side of (3.31) has been broken up into two parts because the 
first part is symmetrical, while the second part, in general, is not except when 
Tr) = Te, and r3 = 7%. 

Let us now consider the expression 


SE[(nyui — --- )*(nw;— --+ )(myuz — «++ (nv; — --- YY. (3.32) 


This is a double summation in which c,;; = c;; and in which the diagonal terms, 
Cii, are Missing. 

Consider next a general term of k factors from the expansion of each bracket 
of (3.32). As we are dealing with symmetric functions, there will be no loss in 
generality if we consider the first k subscripts only; and if we let the lower limits 
of the exponents of the u’s and v’s begin with zero we may consider that each 
parenthesis of a given bracket contributes exactly k factors. Such a term, 
omitting the coefficient, may be written as follows: 


’ ‘ ’ ’ k ‘ 
z 6 > ’ 
E(ug? ++ ugtogt --- oftugt ++. ugwit ++ ote) = [] E(uyettwir*?,) 
h=1 


k 
= |] Peon + ax) (Br + 8). (3.33) 
h=1 
This term occurs in every one of the 3nn, brackets of (3.32), having the same 
numerical coefficient in every one of them, which is 


(a !)? (b!)? 


serve et tana 3.34 
Tla, ! la, ! 118, ! 116, ! ( ) 


To obtain the n; coefficient of (3.33) we break up (3.32) into the following partial 
summations: 


E\(nmyu; — --- )*(nv— --- )>(niu; ae )*(ny0; —---)] = El(mm —--- )* 
(ny. — --- )® (mite — +--+ )*(myve — --- Je) + ee) + El(muen — --- , 
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k 
(nwa — --- )? (mu — --- (mm, — ++ YP) + S B| (ra — +. 
i=l 


(nw; —---)® S (nu — --- )*(nw;— --- » | + S [E{(mu;—---- 


mkt 8,7 =k+1 
(ny; — +++ )(miu; — +++ )*(nw; — +++ 
From this equation we get for the total coefficient in n of the term (3.33) the 
following expression : 


k 


s (ny) trent Orr 8 yy 4 in, S [(—29,) “7% 4 (—m)*s7? a] + Cee = 


h,h’=1 h=1 


l- 


The following restrictions on the a’s and 8’s must be observed. 


a +atees t+a=a B, +Be+--- +8, = 
(a) , , , (b) , , / 
Q, ta,+--- +a, =a B, +Be+--- +8; 


(c) ax +o, +6 +8;, #1. 


From (c) we obtain the upper limit of k, namely: t = a+ b. 
Combining the various above equations we finally obtain: 
n n,b 
(n)2@+») S(USViU%V%) = (a!)2%(b!)?_ S S 
d 


n=l aha »,8r,8 =0 
h h 


k , j k oe 
{ S (amyrrrrrrn ty + S [( ~1,)***"* + (—m)%0"*s] + cy 
h,h'=1 


h=1 


k=1 


Th ’ , 
IT Poanta’ (Br+B,) 


/ 


a (3.35) 
Ila; ! Ila, ! IIG; ! ITs, ! 


Turning to the second part of (3.31) let us consider the expression 
E{(mu; — ---) (nw; — +++) (ru; — ++) (nw; — ---) AP BEAYVB™ 


for a given set of r’s. The term (3.33) may also be considered as a general 
term of this last expression; of course, the exponents of the w’s and v’s will be 
different in this case. In order to evaluate the complete coefficient of a term 
like (3.33) we again write; 
, ri r3 r9 r 
SE((nu; — ---)2(mvi — ++ +)*(riu; — ++ +)8(mwy — + -)° AG BAZ B; 
) 4 / saTip™3,72p"4 
= El(nyu, — --- )*(min — --- (nite — --+ F(a — --- )'A, BA. Bs] 
tn? Tenth 
+- Elojtte — --- )*(mive — --- )*(nyu, — --- (nw, — --- Ae BA, B, 


tos + El(nyu. — -++ )*(miv, — +++ er — +++ BC — +++ dP 
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» r r k rT r 
A, B,°A,2, By.) + S El(mu; — --- )*(my; - --- )7A;'B;’ S 
ian j=k+1 


k 
- )8 (nyo; eerie )®A 7B; - S E\ (nyu; im es )8 (nyo; ee aa ys 


31 
A;*B;* S (nu; — ---)*(nv; — --- )7A5'B;*] + S El(nu;—--- ) 
t=k+1 i,j=k+1 
)6(nw; — --- PATBATB). (3.36) 


It is now quite easy to write down the complete coefficient of a term of the 


form (3.33). The numerical coefficient of this term is the same in every bracket 
of (3.36), and is 


(— 1)S,.(a _ 11) ! (a _— i) ! (b _ r3) ! (b _- rs) ! 
1 (3.37) 








Ta, ! Te, !118,! 08; ! 


The coefficient in m, and Aj! B;3 A‘? B’4 is broken up by (3.36) into the fol- 
lowing four parts: 


. a,t+a, +B, tBar a, 
(—n,) * ™ °" ""'471B'34'2B"4 from the first k(k — 1) brackets. 


n 


k 
A4'1pH73 S A’2B"4 vs S ( yratPn A'1Br3 
Ap, Dy, ARID} = —% Ap, Dh 
A= 


h'=k+1 1 


a,+8,, 


k 
| Qn _ § avis, 


h’=1 


from the next k(n — k) brackets. Similarly 


k ’ P & 
im. 8 (—a) * Ph are | Qa - § AB | , from the next k(n — k). 


h'=1 h=1 


And finally: 


n n n n 
IV. S Aj'!Bi3A"2B"4 = SA;)Bi3SAj2B;t — SAV It Blratew) 
1 1 1 ‘ 


4,j=kt+1 


k n k n k 
— S Aj'By3?A;?Bit— S A;'Bi? S Aj;?Bit— S A;,?B,i S A;,'B;? 


h,h'=1 h=1 h’=1 h=1 h=1 


3 k A 

r r r4 

+ 28 A ;1B;3 S A;?B;# = Qrirs Qrore we Q rire) (r3+r4) = Q,1r3 SA;7B;, 
=] 1 


h= h'=1 

k k k k 

— Qrery SAjBi? — S A;'By3A;?Bi4+2 S Aj;'B,? S A;3?B;}, from the 
1 


h,h'=1 h=1 h'=1 


last c* brackets. 
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The restrictions on the a’s and @’s differ from those given above in that a is 
replaced by a — 7 and a — re, and b by b — r3 and b — 74; and from the restric- 
tion (c) we get for the upper limit of k, in this case, 


ti = atftyt =a + ae m+? SS r3 + r" 






4 


; . Sa ' 
when Sr; is even, or the greatest integer less then 9 when Sr; is odd. 


Combining (3.37) with Cy, --- Co. we get for the general numerical coefficient 
in the expansion of the second part of (3.31), the expression 


(—1)Sr; (at)? (1)? 
Ir,! We, ! Wa; !118,!118;! 
















By an obvious manipulation we have 


k jiu 
I+1+II+IV= 8 | (—n) meen ener Os —1|appeay 44 Oo 
h,h'=1 
k k +s 
S | (my * 1 |ARBI + Qua S | (—my"™ , 1 [aver 
h=1 h=1 


k k Pee k 
— S AtBr § | (—my"™ ‘1 appz S AUB 


h=1 h=1 h=1 
k 


aa! ‘ 
S | (-m)"* ‘ | A;7B;4 + Qrirs Qrorg = Qori+ra)(r34+r4) : (3.38) 


h=1 


Finally, combining the various equations we get the formula: 





n a,b t 


oMp., = N-Poar — pr, + 2(n)-2¢t+» (a!)2(b!)2 S S S 


’ 


ih=1 aha ',8h,8 =0 k=1 
A h 








hyh'=1 h=? 


\ s Se i ee ee ——— 





" a,b Sr . 
ITP (ax, + oy) (Bs + B,) tp a(n )—2(atb+)) (a 1)2 2(b 1)2 Ss (—n)s LSr; 
Ila, ! ITB), ! Ila, ! 11g, ! : hd TysT or TT =O Ir; ! 


»B,7,45 t k + Brt+ a+ 3, 

S S4 S$ [(—m)" " 1) AP BPA™B™ 
| h,h'=1 

aha Bh 8 =O 


k k k 
— § [(—n)t% — 1) APBD S ADB — § [(—m)a,+8, — 1] 


h=1 h=1 h=1 








k k 
ABrt S Aj'Br? + Qrory S [(—mi)otee —1] Aj Bz? 


h=1 h=1 
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k «+p. 
+ Qrirs S [(—m) ii —1) A Hi oo + QrirsQrors —_ Qiri4re) (r3+1r4) 
h=1 


MP (s'+a") (Bat 8) (12) 
Ha,! 118,! Ta, ! 118, ! 


In case the n populations are identical the second part of (12) must vanish, 
and in the first part the summations -° 


‘ 1] i KYCy TE Pranta’) (snt8') 
(ante) (Bits) = 
miei. A!h!---! 





where 1, l2, -- - 1, are the number of repetitions of the pairs of integers 

(a1 + a) (8: + Bi), «++ (ax + a%) (Bx + Bi), respectively. 

We then have the following 

Corollary: The mathematical expectation of the variance, »m,,,, of the product 
moment, pas, in samples of n from a single infinite population is given by 


b,b 


a,a,0, 
2M yar = Brass — Boy + 2n)-2*” (aN)? S 
ha nBrrB =O 


t k! On k ant Bata. +8., ; 
Ss ea Ss (—m)rrr wn 8 [(—my)ante 
aoe etal *** &! 


h,h’=1 h=1 


: 
- II P(e. + a) (Ba + By) 

- at h co 1 : eg 

Pre ss | Te, ! 1B, ! He, ! 118; ! = 





4. The Formula for ,M>,,. Formula (12) can by no means be used mechan- 
ically. It does, however, summarize to a great extent the details in finding 
2M, for any given values a, b. Formulae for 2Mp.,, 2Mp3, have been ob- 
tained, but the one for »M>,, is too long to be included in the paper, especially 
since with a little work it can be easily derived by applying (12). The one for 
2M>r, is given immediately below. 


2Mo, = n-*{nin3S[Pi2 — (P21)*] + n2S[PioPo2 + 4(P3oPi2 — mP31Pi1)] 
— 2nnsSP32P}o + (nz + 2)S(P2oP2oPo2 + 8P20PiiPi1) + 6SP20P2oPo2} 
+ 2n*{nn?S(Pi ,B; + 2Pi,A; — Pi, Pi, B: — 2P3,Pi,4ad 

— 4nnS(Pi,B:Pi, + Pi,A:Pi,) — 2mS[nsP§,B:Piy + 2(2n. — 3)P3},A:Pi,) 
4 6nSPi,Pi,A; + 4nS(Pi,P!,B; + PigPi2A; + Pig A.Pin + 2Pi,Pi,B; 
+ Pi.P},A))} + n{nzS[Pi Bi — (P2.Bi)*] + 4SP3oP2 (Bi + B,) 

+ 3(n5 + m)SP3.A} + 48P2oP)2(Ai + Aj)? — 2nsS[P3 oP 524; 
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4 2P:,Pi,(A; + A)*] + 16SP:,A,;Pi,A; — 4n2S(P!,A,)? 
+ 4(2n} + m)SP;, A,B; — 4n,SP},A;B;P), — 8nsSPuPi, A,B; 

4 8S(P:,B,Pi,A; + P',A,Pi,B,) — 4n2SP‘,A,P:,B; 

— 2nynen'S(QaoP 32 + 2QnuP 31) + 2nenS[6QuP} Pio 

+ Qu(P§oPée + 4P{,P{,)]} +2n-*{nmS(QPi, AB? + 2Pi2A! + 5P$,A2B) 
— mS[Qo(P?,B: + Pi2Ai) + 2Qu8P3.B; + 2P3,A))]} + n-*{n?S[P5. A‘ 

+ 4(P2,A7B? + Pj, A7B)) — 2nS[(Q2A:Bi + QuAz)Piy + QuP3 AB; 

+ QeoP 247] + SIQ2oPi2 + 4Q20(QuPi: + Q11P20)}}. (13)" 








CuapTeR IV. The Mathematical Expectation of the Third Moment of pi 





1. The Mathematical Expectation of 3m»,,. Following the notation of the 
last chapter we shall denote the third moment of pi: about its mean by 3mp,, and 
the mathematical expectation of 3mp,; by 3M». We have then by definition. 


3mp,, = {n—S(az; — z)(yi — y) — pu}', 





and by a well known formula we have: 

Moy, = pi, — &2Moudn — §3,. (4.11) 

The last two terms of (4.11) are given by (1) and (11). To evaluate psy we 

write: 

Bi, = E{nS(x; — 2)(yi — y)}* = n*SE(a; — (ys — y)? 
+ 3n“SE(x; — x)*(yi — y)*(z; — x)(y; — y) 

+ 6n*SE(x; — x)(y: — y)(a; — x)(y; — y) (te — z)(ye — y)- 


—2- 


The first term is simply n—*p33 which is given by (10). The evaluation of the 
second term is not essentially different from the evaluation of the left hand side 
of (3.22), and since all details have been given there we shall omit them here. 
To evaluate the last expression let us write: 


















SE(a; — xz)(y: — y)(z; — x)(y; — y) (ae — x) (ye — y) 

= SE[(U; + A,)(V; + B)(U; + A)(V; + B)(Ui + Ad(Vi + Bi] 

= SE(UV,U;V,;U.V,) + SE(U:V:U;V;U.Br) + --- + SE(A,B;A;B;A,B,) . 
(4.12) 


12In case the n populations are identical this reduces to one of Pepper’s formulae, 
Biometrika, Vol. XXI, p. 238, Cor. 1. 




















rn 
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As there is a great deal of similarity among the various terms of the right hand 
side of (4.12), it will not be necessary to go into the details of the expansion of 
every one of them. We shall, therefore, indicate the details for the expansion 
of only two of them—one symmetrical and one non-symmetrical; and as the 
first two terms are of that type we shall use these for the purpose of illustration. 

Using the u, v notation we have 


SE(U iV U;V;Ui.Vi) = n*SE[(mu;i — --- )(nwi — --- (mu; — --- ) 
(ny; er )(nyux eb eG ) (nyu, ian -++)]. 


The maximum number of subscripts appearing in any term evidently being 3, we 
can write without any loss in generality: 


OO a a ee 
ee ee 
eS a ee. eo 
ily so ym + 8 (rm ee 


(mm. — oo-) + (nw. — -++ )(nyve oa oe) + nus — ~++ )(nivs — --+)} 
S(myu; — --- (mv; — «++ (mu; — «++ (ma; — ++) + SE{(mu; — +--+) --- 


(nm — ---)}. (4.13) 


The coefficients of the various terms arising in this expansion can now be 
found quite easily. For example, the coefficient of P};,, which is, of course, the 
same as the coefficient of P3;, is easily found to be 


2 
ot 4 nttet 4-4 ES 


Nsnns _ NNN 3m — 2) 
- 6 
To evaluate the summation SE(U;V,;U;V;U;Bx) = n~*SE[(mu; — ---) 


(nv; — --- ) (mu; — --- )(mw; — --- )(mw, — --+ )By], we break it up into 
partial summations as follows: 


SE|(nu; — --- )(nv; — --- (nu; — --- (ny; —--- ) (rug — +--+ ) By) 
= E{(nm — --- (ny, — --- )[(myue — --+ )(nive — --+ ) (mus — --- ) Bs 
+ (nie — --- )Bo(nyws — --- )(nws — ---)) + (rim — --- )Bilniee — --- 
(ny. — --- )(nyus — --+ )(mvg — ---)} + Ef{(mim — --- )(mm — ---) 
[(mie — --- )(mwe — ---) + (mus — --- )(nyws — ---)) + (nite — --- 
(nye — --- )(nyug — -++ )(myvg — --- )} Srv; — --- )B + E{(nm — --- 


(nw: — --- )[(me — --- )Be + (nus — --- )Bs] + (nyu, — +> )(nw2 a re 
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(miu: — --- )Bi + (mus — --- ) Bs] + (mius — --- )(niv3 — --- 


ee )B:]} S(mu; dein aa0% 

+ E{(mvi — «++ (mi — +++ )+ (mite — +--+ )(mivz — «++ ) + (mitts — «++ ) 
(nyw3 — --+)} Sra; ee ee ee, ee 
Se en )Bs} Sm i i ae 568 

bite mB ose ES; ee a er 


(nw; — --+ )(myuz — --- ) By. (4.14) 





The expansion of (4.14) is not as difficult as it appears for only two subscripts 
can appear in any term: the explicit appearance of the subscript 3 is due to the 
fact that we are dealing with a triple summation. We, consequently, do not 
need to expand those parentheses in which B appears. 

We shall now, without any further details, state the final result, which is: 


sMpy, = n-*{S[ni P33 — P3oP 33 + 3n,(P3iPi2 + P2oPis) + 3m (ni + 2)P22P iy 

— 3(2nt + 1)P2;Pi, + 3nsPi:P32P 20 + 6(n} + 3n, — 2)P},Pi,P ii) 

— 3n,SP},[S(n{ P22 + P2oPd2 — ni (Pii)? + 2Pi:Piy)] — ni(SPi,)")} 

+ 3n{S[n}(P3.B; + P23A,) + 2a(P2,P},B; + Pi2Pi,A)) 

— InP; Pi LB, + Pi.PiA) — 2n,(Pi PLB, + Pi,PiA) 

+ (Pi2P 2B; + P2:P52A;) — 2n,(Pj.P30B; + P2i:P52A)) 

+ (P5oP 32B; + PosP30A;)}} + 3n—*{Sln,(P3,Bi + Pj;A7) 

+ n(P2oP1,Bi + PosPi:Ai) — (PiP20Bi + PosPi:Ai) 

— 2(P3,B,Pi,B; + P5.A;Pi,A;) + 2n,P:,A;B; — 2P},B;Pj.A; 

— 2P},A;Pj,B; + 2n,P},Pi,A;B; — 2(P},)?A,;B,)} + n{S[(P3,Bi + P3343) 

+ 3(P3,A;Bi + Pj,AiB)]}. (14)" 
Where a = ni +141. 

This formula is shorter and simpler than the formula for 2M>,,, although they 
are of the same order. This is due to the symmetry of 3M»). 
CHAPTER V. Product Moments of Trivariate and Quadrivariate Populations 


1. Some additional definitions and notation. In this chapter we shall indicate 
briefly how the method of the previous chapters may be extended to populations 


13 Cf. Biometrika, Vol. XXI, p. 253, formula (19). 
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of more than two variables. We shall do this by deriving some of the simpler 
formulae, corresponding to those of Chapter II, for trivariate and quadrivariate 
populations. 

The notation will be slightly changed in that we shall symbolize the new 
variables by priming the symbols for the variables used in the previous chapters. 
Thus, we shall indicate the kt" trivariate population by (X;,, Y,, X;) and the 
kt» quadrivariate population by (X,, Y;, X;, Y;), and samples from such 
populations by (2, y¥;, ar,.) and (a,, Yrs Xe Yu) Yespectively. 

We shall denote by P’};, the product moment of the m'® population of order 
tin X,jin Y, and kin X’, and by P?;;,, the similar product moment for a quadri- 
variate population. These are defined by the following equations: 


Tak _ E(X,, -_r On) (Ym prt b,,)(X;, = a)’, (5.11) 
mint = E(Xp — Ap)'(Yin — Bm)! (Xin = Cm) (Ym = dm)! (5.12) 


where am, bm, ete. are defined as in Chapter I part 2. 
The sample product moments corresponding to P?;,, P’,,; will be denoted 
by pijx and pi; respectively. They are defined by: 


Pijk _ n-* S (s,.. cml £)'(Ymn aa y)' (x, on a7’, (5.13) 


m=1 


Piz = n- S (x, - r)*(Ymn ic y)*(a,, —zx’)*(y,, — y’)' . (5.14) 


m=1 


Finally we shall designate E(p ij.) and E(pijm) by pij, and pijm respectively. 


2. The Mathematical Expectation of p,,, and po. By definition we have 
Pu = ElnS(2; — x)(yi — y)(x; — 2’)). (5.21) 
Applying the transformations (1.17) this equation becomes 
npy, = E[(S(U; + A)(V; + B)(U; + C)] = SE(U,V,U;) + SE(U,VC) 
+ SE(U,U;B, + SE(V,U;A,) + vanishing terms + SE(A,B,C,). (5.22) 


Since EA;B:C; = AiBiC;, SE(A:BiC;) = SA;:B,C;. Following the previous 
notation we shall put SA;B,.C; = Qin. 

When the expression SE(U;V;U;) is expanded, no other non-vanishing terms 
except those of the form E(u,v,u;) = Pi,, can appear. The coefficient of this 


term will evidently be the same as that of P}, in (2.23), namely: n-?nyne. 
Whence: 


SE(U,V,U;) = nnn, SPi,. 
The three terms following the first of (5.22) are by (2.24) equal to 
n no S(P 310; + Pio B; + P31: A)). 





